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Absrtract. Group classification of systems of two coupled nonlinear reaction- 
diffusion equation with a general diffusion matrix started in papers pQ, [2] is 
completed in present paper where all non-equivalent equations with triangu- 
lar diffusion matrix are classified. In addition, symmetries of diffusion systems 
with nilpotent diffusion matrix and additional first order derivative terms are 
described. 



1 Introduction 

Group classification can be considered as one of basic elements of qualitative 
analysis of entire classes of differential equations. Such classification makes it 
possible to select the basic non-equivalent equations in the considered class of 
them. In addition, it opens the way to application of powerful group-theoretical 
tools for searching for conservation laws, construction of exact solutions, group 
generation of families of solutions starting with known ones, etc., etc. 

In the present paper we complete the group classification of systems of reaction- 
diffusion equations with general diffusion matrix performed in papers pQ, |2j. 
These equations can be written in the following form: 



u t - A(A u u + A l2 v) = f\u,v), 
v t - A(A 21 u + A 22 v) = f 2 (u,v) 



(1) 



where u and v are function of t, x 1} X2, ■ ■ ■ , x m , A 11 , A 12 , A 21 and A 22 are real 
constants and A is the Laplace operator in R m . 

Equations (pQ) form the grounds for a great many models of mathematical 
physics, biology, chemistry, etc., which makes their group classification to be 
especially relevant. We will not discuss here the history of investigation of sym- 
metries of systems (pQ) which can be found in pQ, [2]. 

Up to linear transformation of dependent variables there exist three a priori 
non-equivalent versions of equations (pQ), corresponding to a diagonal, triangular 
and square diffusion matrix. The equations with diagonal and square diffusion 
matrix have been studied in papers pQ and |2j correspondingly were the complete 
group classification of such equations is carried out. 



Here we consider the remaining case when the diffusion matrix is triangular 
and system ((TJ) is reduced to the following form: 

u t - aAu = f l {u,v), 
v t — Au — aAv = f 2 (u, v) 

where a is a real constant. 

In the case when a = the diffusion matrix is nilpotent and it is naturally to 
generalize (J2J) to the following system 

Ut-ppV Xli = f 1 (u h u 2 ), 
v t - Au = f 2 (u 1 ,u 2 ) 

were (and in the following text) summation from 1 to m is imposed over the 
repeated Latin indices, p^ are arbitrary constants. Using linear transformations 
of independent variables we reduce values of p M to 

Pi = P% = ■ • • = Pm-l =0, p m =P (4) 

where p = a/p? + P2 "I ^ Pm- 

We notice that in addition to possible application in mathematical biology 
equations can serve as a potential equation for the nonlinear D'alembert 
equation. 



(2) 



(3) 



2 Symmetries and determining equations 

To describe symmetries of equations (j2j) and Q with respect to continuous groups 
of transformations we use the Lie infinitesimal approach. Using the standard Lie 
algorithm we can find the determining equations for coefficient functions rj, £ a , ir b 
of generator X of the symmetry group: 

X = V d t + i v d Xv - tt 1 ^ - n 2 d v (5) 

The first set of determining equations describes dependence of rj, £ u and 7r b on u 
and v: 

Vu = vv = o, c = C = o, Ku = < v = C = o. (6) 

So from (jSJ r\ and ^ are functions of t and x M , and 

n a = N al u + N a2 v + B a , a = 1,2 (7) 

where iV afe , B a are functions of £ and x v only. The remaining determining equa- 
tions in the case a ^ are pQ: 

2^ = -5^(rj t A + [A, N]), Vx „ t = 0, (8) 
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Q - 2aNH ~ aA ^ = °> = 0, Nil = ~ aN xli 

Vt fk + N kbjb + ( N u, _ AA ks N sb )u b + B k - AA kc B c (9) 
= (B a + N ab u b )f k a . 

Here N and A are matrices whose elements are N ab and A ab , S ab is the Kronecker 
symbol, and we use the temporary notation u = u\, v = u 2 . 

Using (J?!-® we find the general form of symmetry (jSJ) admitted by equation 

©: 

X = W^x^ + vd t + Pfl d Xft + XK + a^G^ + + fiD , , 

-c\ud u + vd v ) - c 2 ud v - B l d u - B 2 d v 1 1 

where the Greek letters denote arbitrary constants, B 1 , B 2 are functions of t,x, 
and C 1 , C 2 are functions of t, 

K = 2t(td t + ay^J - ^ (\{ud u + vd v ) - ud v ) - tm(ud u + vd v ), 

= td x ^ + \x^ i} a {ud u + vd v ) - jjiud v ) , , , 

G/j, = (d Xll + ±7^ {\{ud u + vd v ) - ±ud v )) , 1 ' 

D = td t + \xpd 9vi . 

For a = symmetry X again has the form (fTUj) where however A = cr M = 
c<J M = 0. In addition, B 2 can depend not on t, x only, but also on u. 

In accordance with (jHJ, (jZj) and © equation (J2J) admits symmetry (fTUj) iff the 
following classifying equations are satisfied: 

(A(m + 4)t + fi + ± {\\x 2 + + 7 e^u; M xJ) + C 1 ) J 1 
+C/n + B, 1 - aAB 1 = (B l d u + B 2 d v + C\ud u + t><%) + C 2 ud v 
+\mt(ud u + + (\\x 2 + + je^u^x^) (~{ud u + vd v ) - jzud v )) f 1 , 

(A(m + 4)i + a* + C 1 ) / 2 + C 2 / 1 + (|Ax 2 + a,x, + ye^u^) (±/ 2 - f 1 ) 
+C}v + C 2 « + B 2 - AB 2 - aAB 1 = (B 1 d u + B 2 d v + C^udu + vd v ) + C 2 ud v 
+Xmt(ud u + vd v ) + {\Xx 2 + a^x^ + ye^u^x^) (\{ud u + vd v ) - ^ud v )) f 2 . 

(12) 

Equations (fT2|l are nothing but a special case of the generic classifying equa- 
tions (2.9), ref. pQ which are valid for arbitrary invertible diffusion matrix. 

Equation (JHJ) needs a particular analysis. For p ^ the related symmetry 
operator (j5J) reduces to the form 

X = X + X (13) 

where 

X = Xd t + v a d Xa + V kl x k d h 
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X = fi {3td t + 2x u d u - vd v ) - F {ud u + vd v ) - B^u - B 2 8 v . (14) 

Here ^>^ u is an antisymmetric tensor and summation is imposed over the 
repeated Greek and Latin indices indices from 1 to m and from 1 to m — 1 
correspondingly. 

Symmetries Xq are admitted by equations (j4*j) with arbitrary non-linearities 
f 1 and / 2 while the classifying equations generated by symmetries X are 

(3 f i + F)f 1 + F t u + B}-pB 2 Xm 

= {B 1 d u + B 2 d v + Fud u + (F + n)vd v ) f\ 

(15) 

(4/i + F)f 2 + F t v + £? 2 — AB 1 

= {B l d u + B 2 d v + Fud u + (F + fj,)vd v ) f 2 

where F and B 1 , B 2 are unknown functions of t and t, x respectively. 

The determining equations for symmetries of equation with p = are 
qualitatively different for the cases, when the number m of spatial variables 
xx,X2, - ■ ■ x m is m — 1, m = 2 and m > 2. The related generator (jSJ) has 
the form 

X = ctD + (f{N — M)dt) d t + 2mH a d Xa - (N + (m - 2)fl£j ud u 

- (M + (m + 2)iZ£j - -B 1 ^ - B 2 d v - B 3 ud v [ } 

where summation from 1 to m is imposed over repeating indices, the Greek letters 
denote arbitrary parameters, M, N are functions of t, B 1 , B 3 are functions of t, x, 
B 3 are functions of t, x, u, and 

H a = 2XbXbX a — x 2 X a if m > 2. 

For m = 2 H a are arbitrary functions satisfying the Caushy-Rieman condi- 
tions: 

rxl tt2 ttI tt2 

XI £2' X2 xi ■ 

In the case when m = 1 if 1 is a function of x and the sums with respect to 
a in (fTo^l are reduced to the only terms. 

The related classifying equations have the form 

(a + 2N — M + (m — 2)H£ a ) f 1 + N t u + B\ = (B 1 d u + B 2 d v 
+B 3 ud v + (N + (m - 2)iz£) w<9„ + (Af + (m + 2)#£j J 1 , 

(a + iV + (m + 2)fT*J / 2 + + M t « + B t 3 « + 5 2 - ABi ^ 
+(2 - m)AH«u = {B 1 d u + S 2 ^ + B 3 ud v + (JV + (m - 2)fT»J «<9 U 
+ (M+(m + 2)//«J^)/ 2 . 

We notice that in this case symmetry classification appears to be rather com- 
plicated and cumbersome. Nevertheless, the classifying equations can be effec- 
tively solved using the approach outlined in the following sections. 
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3 Classification of symmetries 

Following pQ we specify basic, main and extended symmetries for the analyzed 
systems of reaction-diffusion equations. 

Basic symmetries form the kernel of the main symmetry group and so are 
admitted by equation (J2J) with arbitrary non-linearities f l and f 2 . They are 
generated by the following infinitesimal operators: 

P = d t , P x = d Xx , = Xftd^ - x v d Xll (18) 

and correspond to shifts of independent variables and rotations of variables x v . 

Main symmetries form an important subclass of general symmetries pOjl 
which correspond to A = a v = uj v = and so have the following form 

X = HD- C\ud u + vd v ) - C 2 ud v - B 1 d u - B 2 d v . (19) 

To describe all Lie symmetries admitted by equation (J2J we follow the proce- 
dure outlined in [T] which includes the following steps: 

• Finding all main symmetries ()19)) . i.e., solving equations (fT2"j) for = 

v = Pv = Ov = u v = 0: 

(/x + C l )f l + C l tUl + B\ - aAB 1 

= {C\ud u + vd v ) + C 2 ud v + B x d u + B 2 d v ) f\ 

(20) 

(H + C l )f 2 + C 2 u + C\v + B 2 - aAB 2 - AB 1 

= {c\ud u + vd v ) + c 2 ud v + b^,. + B 2 d v ) p. 



Specifying all cases when the main symmetries can be extended, i.e., at 
least one of the following systems is satisfied: 

af 1 = (a(ud u + vd v ) - ud v ) f 1 , , , 

af 2 ~f l = (a(ud u + vd v ) - ud v ) f 2 - 



aif 1 + yu) = (a(ud u + vd v ) - ud v ) f\ 
a(f 2 + yv) - 7W = {a(ud u + vd v ) - ud v ) f 2 

or if equation (j^Tj) is satisfied together with the following condition: 

(m + 4)/ a = m(ud u + vd v )f a , a =1,2 (23) 



If relations (|21jl . (|2*2*|) or (|2Hj) are valid then the system admits symmetry 
G a , G a or K correspondingly. 
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• When classifying equations (J3J) or (J2J) with a = the second step in not 
needed in as much as in accordance with (fTBj) and (fTH ) these equations 
admit basic and main symmetries only 

We note that the first step of the described procedure appears to be very com- 
plicated. In spite of that equations (J2"U|) can be effectively solved using separation 
of variables, the tree of versions of such separations is rather large and includes 
a lot of intersections. 

In the next section we present specific tools used to overcome these difficulties. 

4 Algebras of main symmetries for equation (E3) 

In accordance with the plane outlined in Section 4, to make symmetry classifica- 
tion of equations we first describe the main symmetries generated by operators 
|TT)|) and then indicate extensions of these symmetries. 

First we note that for any f 1 and f 2 equation (j2J) admits the following equiv- 
alence transformations 

u^^u + b 1 , v -» K 1 v + K 2 u + b 2 , 

f^X 2 K'f\ f 2 ^\ 2 (K'f 2 + K 2 f), (24) 
t — > A 2 t, Xb — > A 1 x\ ) 

where K 1 , K 2 and A are constants which are distinct from zero, b 1 and b 2 are ar- 
bitrary constants. In accordance with its definition, equivalence transformations 
keep the general form of equation (J2J) but can change the concrete realization of 
their r.h.s.. For some non-linearities f 1 and f 2 there exist additional equivalence 
transformation which will be specified in the following. 

We will use transformations fl24(l to simplify generators ()19j) . 

To solve rather complicated classifying equations (J2DJ), we use the main al- 
gebraic property of the main symmetries, i.e., the fact that they should form 
a Lie algebra (which we denote by A). In other words, instead of going throw 
all non-equivalent possibilities arising via separation of variables in the classify- 
ing equations we first specify all non-equivalent realizations of algebra A for our 
equations up to arbitrary constants and arbitrary functions. Then we easily solve 
classifying equations (J2UJ) with known functions C k and B a . 

Consider consequently one-, two-, • ■ ■ n-dimensional algebras of operators (|2"5j) 
which we write in the form 

X = fiD + N, N = C\ud u + vd v ) +C 2 ud v + B 1 d u + B 2 d v . (25) 

Let (f2"5j) be a basis element of a one- dimensional algebra A then commutators 
of X with P an d P a are equal to a linear combination of X and operators (jl8j) . 
This can happen in the following cases: 

C a = 0, B a = e xt+u) - x u a , fl = ' 
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where fi ab , fi a , A, and lo = (u>i, u 2 , • • • , uj m ) are constants, to • x = u u x u and uj ^ 0. 

To classify all non-equivalent symmetries ()25|). pUj) we use isomorphism of 
the related iV with 3x3 matrices of the following form 



u 1 fi 1 j . (27) 

V 2 /I 2 fl 1 

Equivalence transformations (j2^j) generate the following transformation for 
matrix (|2j 



g - g' = E^CT 1 (28) 

where 

/ 1 \ / K 1 \ 

[/ = fe 1 , ET 1 = — -6 1 1 . (29) 

\b 2 K 2 K 1 J K \ b l K 2 -b 2 -g 1 / 

Up to transformations (}2*Hj) there exist six non-equivalent matrices g, i.e., 



£i=010, g 2 




g 4 = u i u , # 5 



Let us denote 

$ = ^ 22 m<9 u + g 33 vd v + g 32 ud v + g 21 d u + g 31 d v 
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(30) 



where g ks are elements of matrix g. Then in accordance with (j2*Uj) the related 
symmetries ()25)1 have the form 



X = /j,D + g forg = g 1: g 4 , g 5 , 
X = e xt+ ^~g for g = gi ,g 2 , (31) 
X = e xt g for any g fl3Dj) . 

Formulae (|3T)|) - (pTTf give the principal description of all possible one-dimension 
algebras ^4 which can be admitted by equation (J2J). 

To describe two-dimension algebras A we classify matrices g ()27j) forming 
two-dimension Lie algebras. Up to equivalence transformations tfl^i there exist 
six such algebras: 

^2,1 = {#3, h}, A 2 ,2 = {gugs}, A 2 ,z = {#5,<?2}, A 2,i = {#6, £2}, (32) 
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^2,5 = {01,02}, ^2,13 = {91,93} 



(33) 



where 92 is matrix g 2 (j3*Uj) with A = 0. 

Algebras (}3*2j) are Abelian while (|33J) are characterized by the following com- 
mutation relations: 



[ei, e 2 ] = e 2 . 



(34) 



Two-dimension algebras A generated by (J32|) and (|3*3*|) are spanned on the fol- 
lowing basis elements 



and 



< \iD + it + vte 2 , e 2 >, < \lD + e 2 + istei, e\ > 

< (j,D - ei, uD - e 2 >, < F1&1 + G\e 2 , F 2 e,\ + G 2 e 2 > 



< fj,D - ei, e 2 >, < /iD + ei + z4e 2 , e 2 > 



(35) 



(36) 



respectively, where {Fi,G±} and {F 2 ,G 2 } are fundamental solutions of the fol- 
lowing system 



F t — XF + aG, G t = aF + 7 G 



(37) 



with arbitrary parameters A, a, a, 7. Arbitrary parameters \x and v in particular 
can be equal to zero. 

In addition, there exist two dimension algebras A which are induced by one- 
dimension algebras of matrices 9 (J2"?J) , namely 

(38) 



<Fg,Gg>, <fiD + \e vt+ "- x g, e vt ^ x g > 



with F and G satisfying (J37|) . Such algebras correspond to incompatible classi- 
fying equations (J2U|) . 

Up to transformations (J2"8j) there exist four three-dimension algebras v4 3jl — 
A 3i 4 of matrices ()27|) and the only four- dimension algebra of such matrices which 
we denote as v4 4 : 

Table 1. Three- and four-dimension algebras of matrices (1271) 



Algebra 


Basis elements 


Nonzero commutators 


^3,1 
A 3 ,2 
A 3 ,3 

A, 


ei = g%, e 2 = 03, e 3 = g 2 
ei = 05, e 2 = gi, e 3 = g 2 
ei = 02, e 2 = g 5 , e 3 = g 6 
ei = 02, e 2 = gs, e 3 = # 4 

ei = 0i, e 2 = # 3 , e 3 = # 2 , 

e 4 = 05 


[ei, e 2 ] = e 2 , [ei,e 3 ] = e 3 
[e 2 ,e 3 ] = e 3 
[e 2 ,e 3 ] = e x 

[ei,e 2 ] = e 2 , [ei,e 3 ] = e 2 + e 3 
[ei,e 2 ] = e 2 , [ei,e 3 ] = e 3 , 
[e 4 , e 2 ] = e 3 
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Using commutation relations present in the table we come to the following 
related three-dimension algebras A : 

< \iD - ei, e 2 , e 3 >, < ex, F x e 2 + Gie 3 , F 2 e 2 + G 2 e 3 > 
with e a belonging to A 31 ; 

< /iD — 2ei, vD — 2e 2 , e 3 > 
with e a belonging to A 32 ; 

< fiD — 2e 2 , z/D — 2e 3 , ei >, < ex, D + 2e Q + Ivtex, e a > >, 

< e^+^ei, e^' x e a , e Q / > 

where a, a' = 2, 3, a' ^ a and e a belong to ^4 3i3 ; 

< fiD - Tex, e 2 , e 3 >, < e x , e ut+ ^ x e 2 , e vt+ ^ x e 3 > 

with e a belonging to t4 3 4 . 

The four-dimensional algebra A 4 induces algebras A given below: 

<fiD- 2ex, vD - 2e 4 , e 2 , e 3 >, < e^^ei, e ut+UJ - x e 4 , e 2 , e 3 > . 

Thus we had specified algebras of main symmetries which can be admitted 
by equation (|2j). 

5 Solution of classifying equations for the case 
of invertible diffusion matrix 

Applying results of the previous section we can easily classify main symmetries 
of equation (J2J). Such classification reduces to solving equations ((20} with their 
known coefficients C l ,C 2 and B X ,B 2 which can be found comparing (fT9*j) with 
the found realizations of algebras A. To complete the group classification of 
equations (J2J) we will specify all cases when relations (|2"T|l - ([22jl are satisfied, i.e., 
when the main symmetries can be extended. 

Solving of equations fl20D with known C l ,C 2 and B x ,B 2 is a rather routine 
procedure. We restrict ourselves to presentation of an example of such solution 
in the following. Note that asking for invariance of (J2J) w.r.t. one-dimension 
algebra A we fix the non-linearities f 1 and f 2 up to arbitrary functions while an 
invariance w.r.t. a two-dimension algebra A usually fixes these non-linearities up 
to arbitrary parameters. 

We will solve classifying equations (j^Ujl up to equivalence transformations 
U -> U = G{U,t,x), t -> t = T(U,t,x), x -> x = X(U,t,x) and / -> / = 
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F(U, t, x, f) which keep the general form of equations (J2J) but can change functions 
f 1 and f 2 . 

The group of equivalence transformations for equation (J2J) can be found using 
the classical Lie approach and treating f 1 and f 2 as additional dependent vari- 
ables. In accordance with their definition, equivalence transformations include 
all symmetry transformations (i.e., transformations generated by operators (JTHj) 
and other symmetries which will be found in the following) and also transforma- 
tions (|24j). In addition, for some particular non-linearities f l and f 2 there exist 
additional equivalence transformations, whose list is given in formulae (|H9j) : 

1. it — > exp(ut)u, v — > exp(u)t)v, 

2. u — > u + ut + fix 2 , v — > v, 

3. u — > u, v — > v + pt + px 2 , 

4. u—>u + pt, v — > v exp(pt), 

5. it — > u, v — > v + ptu, 

6. it — > exp(to>t)it, v^v + Ktu + p^-, (39) 

7. m — > it, f — > v — ptu + pAy, 

8. u — > exp(pt)it, t> — > exp(pt) + e^uj , 

9. ix — > u + pt, v — > v + ptu + py, 

10. u — > exp(u;t)ii, f — > exp(atf)(t> — wtit), 

11. Transformations G48I) valid for a = only. 

where $(it) is an arbitrary function of w, iq and F 2 are functions of it which 
appear in the classified equations. 

In the following we specify additional equivalence transformations (|39*|l ad- 
mitted by some of equations (J2J. 

Let us present an example of solving of classifying equations. Consider the 
first of algebras A given by relation (J35j) with ei, e2 belonging to algebra A 2i2 
(|32*J) . It includes two basis elements 

Xl = p-D — M<9 U — vd v , X 2 = vD — ud v . 

Comparing X\ with X (|25jl we conclude that in this case C 1 = 1,C 2 = B 1 = 
B 2 = and so the classifying equations f!2()l) are reduced to the following ones: 

(p+i)r = ( w a u + ^)r, a = 1,2 

General solutions of this system have the form 

fi = u^Fx, f 2 = v^ +1 F 2 (40) 

where iq and F 2 are arbitrary functions of -. 

Thus equation (j2J admits symmetry X\ iff the related non-linearities are of 
the form (pEUj). 
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Asking for symmetry of equation (J2J) w.r.t. transformations generated by X 2 
we come to the following classifying equations (|2U|): 

vf = ud v f\ uf + f^udj 2 . (41) 

Substituting (|4T?j) into (}4~Tj) we come to the equation whose general solution 
for \x 7^ is 

/i = \ u » +1 e yV - , f 2 = e< (Xv + au)u». (42) 

In the special case \i = the solution has the form 

f = \ ue v l + uu, f 2 = e v u (Xv + au) + uov. (43) 

However, in this case there exist the additional equivalence transformation (f3T)|) 
given in Item 1 with p = u, which reduces parameter u in (f4*3*j) to zero. So 
without loss of generality we can restrict ourselves to solutions (}4*2"j) for any p. 

We see that equation (j2J) admits the two-dimension algebra of main symme- 
tries spanned on Xi,X 2 provided f 1 and f 2 have the form ()42|) . This symmetry 
can be extended if functions (jlSJ) and f 2 satisfy one of conditions (|2"T|) . (j22J) or 
both the conditions (jHD , ((2H1) - 

Equation (}2"2"j) is incompatible with (jUJ). In order equation (}2"T]) be satisfied 
we have to impose the condition /i = —au on parameters p, v and a. The related 
equation (J2J), (|4^j) has the form 

v t - Au - aAv = e u ™ (Xv + ou)u~ av 1 ' 



and admits the Galilei generators G M (JTTJ). Finally, asking for equation (j23jl be 
satisfied we obtain one more condition v = — -m which guaranties invariance of 

a & 

equation (|4^|) w.r.t. the conformal generator K of (jllj) . 

The obtained classification results are presented in Table 2, Item 3. 

In analogous way we solve classifying equations for other algebras A and 
specify the cases when the main symmetry can be extended. 

6 Classification results for equations (121) with 
invertible diffusion matrix 

The classification results are given in the following Tables 1-5 when we also indi- 
cate the additional equivalence transformations (AET) (jSHJ) which are admitted 
by some particular equations (J2J). The symbols D, G u , G u and K are used to 
denote operators (JTTJ). In addition, we denote 

K = K + —— (t (pud u + (2 - X)vd v ) + ud v ) . 
A — 1 
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To save a room we present in Table 2, Items 8-11 and Table 3, Items 1-3, the 
classification results which are valid for equation (J2J) with a^O and a = as well. 
The completed analysis of symmetries and the corresponding non-linearities for 
equation (J2J) with a = is given in Section 7. 

In the following tables F\, F 2 and F are arbitrary functions whose arguments 
are specified in the third column, *f>(x) is an arbitrary function of x±, x 2 , ■ ■ ■ , x m 
and i(} u (x) is a solution of the linear heat equation (d t — A)ip u = uip u . In addition, 
we denote by ^^(x) a solution of the Laplace equation A\E f A1 (x) = ^^(x). 

Greek letters in the tables denote arbitrary parameters which can take any 
(including zero) real values. The only exception is parameter e in as much as 
without loss of generality we can restrict ourselves to its values e = ±1. 

Table 2. Non-linearities with arbitrary functions and symmetries for 

equations ((5]) 



No 


Nonlinear terms 


Argu- 
ments 

of 
Fi, F 2 


Symmetries 


1. 


/ = vu + F x , 

f 2 = vu 2 + F lU + F 2 , 


2v - u 2 


ijj v (ud v + d u ) 


2. 


£l ^VU TP 

J = e ri, 

p = e™{F 2 + F lU ), 


2v - u 2 


vD — ud v — d u 


3. 


f 1 = F u 

P = F 2 + uv, 


u 


ip v d v 


4. 


P = au + /i, 

f 2 = vv + F, an = 


u 


e (u-a)t ( u _ ^ Q v 


5. 


f 1 = u\ 

P = uv + uv + F 


u 


e vt ud v , 

e vt (d v + tud v ) 


6. 


f l = (u 2 - 1) , 

p = (u + v ) v + F 


u 


e {»+i)t (ud v + d v ) , 
e^" 1 )* {ud v - d v ) 


7. 


f l = (u 2 + l), 
f 2 = (u + u)v + F 


u 


e vt (cos tud v — smtd v ) , 
e ut (smtud v + costd v ) 


8*. 


p = e uv F u p = e vv F 2 


u 


vD — d v 


9*. 


p = e vu F u 
p = e vu F 2 


V 


uD - d u 


10*. 


P = uu + F u 
f 2 = —afiu + F 2 


V 


e^ +a ^ a (x)d u 


11*. 


p = + vlnu), 

p = v(F 2 + u\nu), v^O 


u 

V 


e ut (ud u + vd v ) 



The items marked by asterisks are valid for both cases a ^ and a = 0. If 
a = then in Item 8* without loss of generality F\ — 1. 
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Table 3. Non-linearities with arbitrary functions, symmetries and 

AET for equations (0) 



No 


Nonlinear terms 


Argu- 
ments 

of 
F\, F 2 


Symmetries 


AET 

(EH 


1*. 


p = uF\ - vv, 
p = ul(v-u) 
+uF 2 -vF 1 , u^O 




e vt (ud v — ud u 

-V8 V ) 

& G a if a = 1 




2*. 


p = u u+1 F h 

p = u u (F 2 u - F x v) 


V 


T / / ) _L r, 1 n 1 s~) 

IS LJ ^ UbUy UjUu 

—vd v & G a 
if v = 0, a = 1 




3*. 


J — a r l? 

p = U » +1 F 2 


V 

u 


fiD — ud u — vdy 


1 if 

fj, = 


4. 


p = e u uF lU , 

p = e < (F lV + F 2 ) 


U 


vD — ud v 


5 if 

v = 


5. 


f l = u(F 1 -u), 
p = F lV + F 2 , 


U 


e vt ud v 

& Tpydy if Fl — V 


3, if 
Fi = 



Table 4. Non-linearities with arbitrary parameters and extendible 
symmetries for equations ([2]) with a^O 



No 


Nonlinear terms 


Main 

symmetries 


Additional 
symmetries 


AET 

(EH) 


I. 


J 1 = An, 
P = au» 


1p0dy, 

e~ xt ud v 


e xt (ud v + Xd u ) 
if n = 2 


3; 

5 if 

A = 


2. 


P = Xe u , 
P = ae u 


D-d u , 


ud v if A = 


3; 

5 if 

A = 


3. 


p = Xu u+1 e^, 

P = e^(Xv + au)u u 


\iD — ud v , 
uD — ud u 

-Vdy 


G a if v — a/i 
& K if v = ± 

m 


1, if 
v = 0; 

5, if 
// = 


4. 


p = \ u »+\ 
p = au^ l+1 , 
Act = 


\iD — ud u 

-Vdy, 
ljj dy 


ud v if A = 


3; 

5 if 

A = 
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Table 5. Non-linearities with arbitrary parameters and 
non-extendible symmetries for equations (J2j) with a^O 





]\[/"i"n 1 1 n t* f*fST*Tn a 

1 N UllllllCdl LCI lllO 




AET 

432) 


1. 


/! = \v»+\ 

J — V>v 


i/D — ud u — vd v , 


2 


2. 


p = Xe ^ p = ae" 


D-d v , V (x)d u 


2 


3. 


p = Ae M , f 2 = aue u 


D-d u - 2ud v , tp d v 


o. 

-J; 

5 if 

A = 


4. 


p = ue { ), 

f = {vu + ii) e A ( W ) 


XD — d v , d u + ud v 


9 


K 
O. 


f 1 = fihiv, 
P = v In v 


^ (x)d u , D + ud u + vd v 
+ ((// - ua)t - ^mi 2 ) <9 n 


o 
z 


6. 


r = x 

f 2 = £ In u, 


D + tt<9 u + vd v +etd v , 
ibnd„, (u — Xt) cL 


3, 7; 
& 5 if 
A = 


7. 


p = Xu v+1 , 
p = \u» +1 Inu, 


vD — ud u — vd v — ud v , 
ip d v 


3 


8. 


Z 1 = e(2w-w 2 ), 

/ 2 = (// + eu) (2u - u 2 ) 

-e^-u, u^O 


X x = ei* (2d u + 2ud v 
+ued v ) , 2tX x + ee^d v 




9. 


f 1 =e(2v-u 2 ) 1 

f 2 = (fx + eu) (eu — u 2 ) 

+E 1 ~fu 


x ± = ei x±\ (2d u + 2ud v 
+e(a±l)d v ) 


9, if 
li 2 = 1 


10. 


f 1 = e{2v-u 2 ), 

f 2 = (fi + eu) (2v - u 2 ) 


e' xt (2cost(d u + ud v ) 
+e(2/icost — smt)d v ) , 
e^(2siat(d u + ud v ) 
+e(2/Li sint + cost)d v ) 





In the following table we present symmetries of a special subclass of equa- 
tions (J2J) whose r.h.s. is given in the table title. The related non-linearities had 
appeared in the classification procedure of reaction-diffusion equations with the 
unit diffusion matrix [2] where they corresponded to a very rich spectrum of 
symmetries. 

In the case of triangular diffusion matrices studied in the present paper the 
classification results present in Table 5 include thirteen non-equivalent types of 
equations. Among them there are five equations invariant w.r.t. two dimensional 
algebras of main symmetries, and four equations admitting symmetries G v i.e., 
being invariant w.r.t. the Galilei group. 
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In Table 5 the following notations are used: 5 = |(/U — v) 2 + Act, u q = 
+ z/), u)± = uJq ± 1 where /i, z/, A and a are parameter used in definition of 
the related non-linearities / , f 2 . 



Table 6. Symmetries of equations with non-linearities 

f 1 = Xv + fiulnu, f 2 = A^ + (au + nv) Inn + uv and a^O 



No 


Conditions 
for coefficients 


Main symmetries 


Additional 
symmetries 


1 


A = 0, 


^ CcJ v^* '^J - 


t/vdu if = 0, 
& G a if ^ = acr 


G a ,if fj, ^ 0, 
acr = ^ — /i, 


2 


a = 0, 
//A ^ 0, 


e"* fAi?9 R + (a- iy)ud„) 


G a if i/ = 0, 
a// = —A 


G a , if a{y — y) = X 


3 


5 = 


+(// - //)«(%) , 
2e u>ot ud v + tX 4 


G a , if u = 0, 
av = -A ^ 


G a , if cj 7^ 0, 
2A = a{n - v) ^ 


4 


A^O, 
5= 1 


(XRd R + (w+ - n^udy) , 
e"- 1 (XRd R + (cu_ - 


G a , if a/x = A 

c<j + c<j_ = 


G Q , if ^+ ^ 0, 
A = — a(uj+ — /i) 


5 


A^O, 
5= -1 


e^AcosiRc^ 

+ ((// — cost — 2 sm£)Mc%], 

e" ot [2XsmtRd R 

— /i) sint + 2 cos£)mc%] 


none 



IfA = /i = 0orA = z/ = then the related equation (J2J) admits additional 
equivalence transformations 10 or 5 from the list (|3*^j) correspondingly. 

We see that there exist a number of non-equivalent systems (0) with non- 
degenerate diffusion matrix. Rather surprisingly the number of such equation 
with nilpotent diffusion matrix (which are classified in the following section) 
appears to be even more large. 

7 Group classification of reaction-diffusion equa- 
tions with nilpotent diffusion matrix 

7.1 Equations with first derivatives in x 

Consider now equations ©, (JH) and specify their Lie symmetries. In this sub- 
section we restrict ourselves to the case p ^ when generators of admitted Lie 
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group have the general form (jl3j) while the related classifying equations are given 
by formula (|15|). Moreover without loss of generality we put p = 1. 

We solve the classifying equations using the technique developed in Sections 
5 and 6. The general analysis of admissible algebras A can be carried out in 
complete analogy with Section 4. Moreover, the results present in Section 4 can 
be extended to the case of equations Q, © provided we make a formal change 
D — > D = 3td t + 2x v d v — vd v in all formulae where the operator D appeared, 
and exclude all algebras A where matrices g^gs and g 6 (1HU|) appear. Of course 
it is necessary to take into account that in contrast with D operator D does not 
commute with d v . As a result we come to the following one-dimension algebras 

X^p = fiD — ud u — vd v , 

x[ 2) = D-vd u , = e» t (ud u + vd v ), 

X[ 3) = D + ud u + vd v + ud v , X® = e a ^ x (d u + d v ), 

X 3 (1) = e ait+ P^ x d u , Xf ) = e^+P^dv 

and two-dimension algebras 

h =< D,X® >, A 2 =< X? ,xf >, A, =< lfUa (1) >, 

A 4 =< X^,^ >, A e =< D + A(ud u + vd v ) + td v , Xf ] >, (46) 

A 5 =< X[ 1] , Xi 2) >, A 7 =< D + 3(ud u + vd^) + td u , x£ ] > . 

In this way the problem of group classification of the equations with the first 
order derivative terms reduces to solving the classifying equations (fT3j) which 
their known coefficients B l ,B 2 ,F and specifying the case when these equations 
have non-trivial solutions. These coefficients are easily identified comparing ()45|) . 
(I46|) with (|T4^) . For example, for symmetry X^ we have F = 1, B 1 = B 2 = 0, for 
Xi the values of these coefficients are B 1 = B 2 = F = 0, fi = 1, etc. Solving 
the related classifying equations (|15jl we easily find the related non-linearities J 1 , 
f 2 which are given in Table 7. 

In accordance with the results present in Table 7 equations 0, with p ^ 
can admit neither Galilei nor conformal symmetry transformations. This result 
follows directly from formulae (fT3j) . (fT4^1. equations. 

In six cases enumerated in the table the corresponding equations ©, (jlj) 
admit infinite dimension symmetry algebras whose generators are defined up to 
arbitrary functions or arbitrary solutions of linear equations, see Items 5-7, 9-14 
here. However, these infinite symmetries generate the same classifying equations 
as one dimensional algebras X\ — X^ and two dimensional algebras A 3 , A A , A 7 . 

Finally we note that in addition to their symmetries, equations (JHI), Q with 
p ^ admit equivalence transformations (J2H) with K 2 = 0. These transfor- 
mations which change functions f , f 2 and also the value of parameter p where 
used to simplify generators (|43|) . (JH3) and the corresponding equations Q, ©• 
Equations 0, © with the non-linearities given in Items 3 (when v = 0) and 8 
of Table 6 admit additional equivalence transformation u — > e at u, v — > e at v. 
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Table 7. Non-linearities and symmetries for equations (JU) with 

p — 1 





Non-linearities 


Arguments 

n f El El 


Symmetries 


1. 


J — u i 

J — U ^P 2 


vu"^ 1 


2fiD — ud u — vd v 


2. 


f 1 = v 3 F u 

f2 _ ,.4p 


u — v In v 


D - vd u 


3. 


J 1 = u(Fi + vlnu), 

f 2 _ rp _i_ I _ A 


V 

u 


e vt (ud u + vd v ) 


4. 


f 1 = ufF u 

/ = «l ^2 


v — v In u 


D + ud u + vd v + vd v 


5. 


Z 1 = Aw + -Fi, 

f 2 — — //7/ -1- 


V 


e xt ^ u (x)d u 


6. 


f l = vv + F 1 , 
f 2 = Xv + F 2 


u 


e xt ~ UXm ^/(x)d v 


7. 


f l = Xu + F x 
p = av + F 2 


u — V 


e e 2 ^fi{x, x m -f- 1) {u u -f- o v ) , 


Q 
O. 


T" f /l 1" 

J — cra-L u 2 , 

/ 2 = ™rM 




D, ud u + vd v 


Q 


f 1 _ ^^311 

/ — ae , 
f 2 = ve 4u 




D - d u , V(x)d v 


i n 


f 1 _ ryp-Zv 
J — CtC , 

f 2 = ue~ 3v 




D + ud u + vd v + d v , ^o(x)d u 


11. 


p = au^ + \ 
p = vu ^+ x 




[iD — ud u — vd v , ty(x)d v 


12. 


p = av 2u+ \ 
p = uv 3u+1 




vD - ud u ~ vd v , ^ (x)d u 


13. 


f L = ctu*, 
P = vhiu 




\l) + vd v + ud u + vtd v , ty(x)d v 


14. 


f l = i/lnv, 

f2 a 




\D + ud u + vd v + vtdu, ^o(x)d u 
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Here ^(x) and ^{x,x m + t) are arbitrary solutions of the Laplace equation 
= ^ ^ in m-dimensional space, ^(x) is a solution of the Laplace equation 
in m — 1-dimensional space, x = (xi,x%,- ■■ ,x m _i), /i, v and A are arbitrary 
parameters satisfying v\ ^ 0. Finally, we denote D = 3td t + 2x v d v — vd v . 

7.2 Equations (0) with a = 

The procedure of classification of equations ©, (jH) with p = (or equations (j2J) 
with a = 0) appears to be more complicated then in the case of p non-zero. The 
general form of symmetry admitted by this equation is given by equation (jTHj) 
while the classifying equations take the form (JT7J). 

A specific property of symmetries (fTBj) is that in contrast with operators (fTD|) 
and (jlHj) they can generate non-linear transformations for dependent variables 
since the coefficient B 3 can be a function of u. In addition, the classifying equa- 
tions critically depend on the number m of independent variables x v and are 
qualitatively different for the cases m — l,m — 2 and m > 2. Nevertheless, 
these equations can be effectively solved with using the analysis of low dimension 
symmetry algebras present in Section 4. 

One more specific point in the classification of equations (J2J) with a = is 
that they admit powerful equivalence relations 

u^u, v — > v + $(-u) (47) 

and 

u — > u, v — >• v + $(u,t,x) (48) 

which did not appear in our analysis presented in the previous sections. 

Equivalence transformation (|4*Tj) (with $(u) being an arbitrary function of u) 
are admitted by any equation (j2J) with a = 0. Transformations (j4~Sj) are valid 
for the cases when f l does not depend on v and in the same time f 2 is either 
linear in v or does not depend on this variable. Moreover, the related functions 
$(w, t, x) should satisfy the following system of equations 

ffit - - f%v = 0, 

Thus the group classification of equation (j2J) with a = is reduced to solving 
the classifying equations (fTTj) with using the algorithm presented in Section 3. 
We will not reproduce here cumbersome and rather routine calculations which 
are needed to classify equations (J2I) with a = which can be carried out in the 
same way as in Sections 2-6 and taking into account the specific points mentioned 
in the above. The classification results are presented below in Tables 8-10 and 
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(49) 



also in Table 2, Items 8-11 and Table 3, Items 1-3. The related items are marked 
by asterisks. The additional equivalence transformations are specified below the 
Tables 8,9 and 10. 

In Tables 8-10 the symbol W denotes a function of t, x and u which solve the 
following equation: 

f v -w t - w u f = 0. 



Table 8. Non-linearities with arbitrary functions for equations (J2J) 
with nilpotent diffusion matrix 



No 


Nonlinear terms 


Argu- 
ments 
of F a 


Symmetries 


1. 


p = F^-", f 1 = F 2 u» 


V 


Q x = {p, — l)D - vtd t - ud u 
-(v + 1)^ & (to - 2)x 2 d Xa 
-x a Qi if v(m - 2) = 4, 
\i(m — 2) = to + 2, m ^ 2 


2. 


p = F lU v»-\ p = F 2 v fM , 
F 2 ^0 


U 


uD — td t — vd v & 
e w d v if fi = 1 k 
H a d Xa - H° h vd v if to = 2 


3. 


P = F lV -\ p = F 2 + uv 


u 


e 1 ^ (9 t + uvd v ) & e*^ if 
F 1 = 


4. 


p~l T — 1 it "1 C 1 ) T — 1 (/ 

f l = Fiv^ , f — F 2 v^ 


ve u 


//£> - td t - vd v + d u 


5. 


p~\ fr\ r ,r > t—\ 

/ - 17 + ^ P = F 2 + vv 


ve u 


Jjf / O O O \ 

e vt (<% + - ^<9 U ) 


6. 


P = 0, f 2 = F 2 


V 


#o(z)<9«, aJ Se„ + 2u du 


7. 


P = F u p = 


u 


e w d v , x a d Xa -2vd v 


8. 


p = -^- i u + F 1 u 2 -^ 
p=^v + F 2 u, n±\ 




((1 - fi)tdt - vud u 
—ufivd v ) 


9. 


f 1 = uFi, m = 1 
P = vF 2 + u 


vv? 


Q 2 = cos(2x) (ud u — 3vd v ) 
+ sin(2x)xd x , Q 3 = (Q 2 ) x 


10. 


f 1 = uFi, m = 1 
p = v F 2 -u 


vu 3 


Q 4 = e 2x (d x + ud u - 3vd v ) , 
Q 5 = e~ 2x (d x - ud u + 3vd v ) 


11. 


p = F h p = vF 2 , m = 2 


ve u 


H a d Xa - H b Xh (vd v - d u ) 


12. 


P = z/e« , p = e«F 


u 


D — ud v 


13. 


P = F U F 2 = vF 2 + F 3 


u 


e w d v 



For the non-linearities enumerated in Items 2 (when u = 1), 3 (when F\ = 
0), 4 and 8 of Table 8 the related equation (J2J with a = admits additional 
equivalence transformations (|48jl. In addition, transformations (|24|) and (H7I) 
and some equivalence transformations from the list (}39|) are admissible, namely, 
transformations (2) for the non-linearities given in Item 1 (when v = — 1, /j, = 0) 
and in Item 6, transformations (1) and (3) for the non-linearities from Item I 
(when v — 1, \i = 0) and Item 7 respectively. 
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Table 9. Non-linearities with arbitrary parameters and extendible 
symmetries for equations (J2j) with nilpotent diffusion matrix 



No 


Non-linearities 


Main 

symmetries 


Additional 
symmetries 


AET 


1. 


p = Xu u+1 v», 


(fj, + u)td t 
-(//+ l)ud u 
+ (u - l)vd v 

Qq = 2fiud u 
+(fj,+ u)x a d Xa 
—2vvd v 


x a Q 6 - 2nx z d Xa if 
k{vtl + 2) = v, 
k{2 — m) = \i 




e w d v if A = 0, 
/i — — 1 & IX o Xa 
-(m - 2)x a Q§ if 

, , m+2 m J o 

V ~ m-2> 771 T * 


J-J-, 
3,5 


^o(x)d u , if 

a = 0,v = -1, & 

X a Q% "1" m -\-2 X ^ x a 

if /u = m ;2 m ^ 2 

" m+2 ' ' 


11, 
2 


if A = v = 


2 


e u v ii fi — u &i 
H a d Xa - H«vd v 
if m = 2 


1 1 J?r 
11 

5 if 

\ = a 


2. 


p = \ u p+l v-\ 
P = ou v + ev, 

A^O 


e et (d t + evd v ) , 

Q = Qq\ii=~1 


X a Q Q m -2 X 

•f . , m+2 / o 

11 V = m-2' m T 2 




^oO^R if 

(7 = U, 1/ = — i 


2 


u d v + ud u 

if cr = 0, i/ = 1 


1 


3. 


j = Ae , 
p = ae^ +l >, 
Xa = 


[y + \)u — vd v 
-td t - d u , e w d v 


i/ [va v + ta t ) - d u 
if a = 


1 1 

11, 
3 


uc^ + £<9 t 

11 A = U 


11, 

d, 


4. 


p = Ae^ +1 >, 
/ 2 = ae vv 


(v — 1)D — ud u 
+td t - d v , 
^o{x)o u 


v (ud u - ^) 
+a„ if A = 


2 

for 

any A 


5. 


p = \ v ^ 1 e u , 
p = av^e u 


D-d u , 

td t + vd v - fid u 


e w if fi = lk 
H a d Xa - J3£ 

IT TY1 V 

11 / / 1/ 


11 if 

H = 1 


6. 


P = Xlnv, 

„n M+l 

f z = av 2 


fJiD - eptdt 
-K^ud u - vd v 
-Xtd u , ^ (x)d u 


u9 n + 2va„ + td t 
+2Xtd u if a = 


2 


7. 


J 1 = av 1 ^, 
P = Xlnv 


\iD + fiud u — vd v 
-(l + fj)tdt 


£> + ud u — td t 
if a = 


2 
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Table 10. Non-linearities with arbitrary parameters and non 
extendible symmetries for equations (J2j) with a = 



No 


Non-linearities 


Condi- 
tions 


Symmetries 


AET 


I. 


fl = Avx^+V, 

J — U LI V Cx LI j 


fi^O, 
m = 1, 
a = — 1 


-(/j + l)t/<9„ 

Q21 Q3 








ft 7= u, 
m = 1, 
a = 1 


Qa, Q5, Q7 




2. 


p = \u~ 2 v-\ 


a = — 1 


pSt f Pi 1 ci^Fi \ 

e ^ct -|- tuu v ) , 
Q2, Q3 


11 if 




P = au~ 3 + ev — au 


7T7 — 1 

a = 1 


c ^c/^ 1 c.uu v j , 
Oa Qc 


A = 


3. 


p = av^- v+1 




(/i — 2v) D + t4d 4 
— ?;(9,, — (v -\- 1 1?//?. 


2 


4. 


P = Xv, p = e- v 


A ^ 


2D - td t + ud u 
+<9,, + Xtd, ^n(x)d. 


2 


5. 


P = Xe\ f 2 = ae v 


Aa ^ 


D-d v , M%)du 


2 


6. 


p = Xu u+1 e^, 

f 2 — e^liXv + ov\u v 

J C- <* I /\ U \^ \J Ui J Uj 


/iA ^ 


fiD — ud v , 
u V) — iih — 

IS l -S Ub U>n U \Jy 




7. 


P = filnv, 
P = ulnv 


i/ ^ 


^o(x)d u , 

D + ud u + vd v 

* 2m r ) Uu 


2 


8 


f 1 — f 2 — fin?/ 


f hi 

C —I L 


D — td t + ud u 
4-rtf) tf), 4- 1)3 

TCt-Cu, v^t \ UKJv, 

<3>(tt, x)d v 


3,5, 
11 


9. 


f 1 = e (In v — k In u) u, 
P = e (In v — k In u) v 


m ^ 2, 

' m— 2 


( 1 K\ r T Ft 

+2nvd v + 2ud Ul 

e (l- K )et ( u Q u + v Q v j 


1 

if 

AC = 1 


10. 


~f ^ — C'? / ^ ^ <m 1 In ii 
J — C(X i 1 //i T Zl 111 [X 

+ (2 - m) Inv) , 


rn =L 1 9 
// 1 y= 1 , ^ 

a = 


e 4et + vd v ) 






/ 2 = £v ((m + 2) In it 

_|_ i v id i In 1 1 i /t/h 

I ^ / / 1 J in t/ j (jf.Lt 


m = 2, 


H«d Xa -Hlvd v , 

p^.St ( q , pi _|_ nff) \ 








m — 1, 
a = 1, 
f = 1 


Q2, Q3, 

e 4 * + vd v ) 








m = 1, 
a = 1, 
f = -1 


Qa, Q5, 

e~ u (ud u + vd v ) 
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Table 10. Continued 



No 


Non-linearities 


Conditions 


Symmetries 


AET 
ii 


11. 


f 1 = fxulnu, 

/ = \iv In u + vv 




e w d v , 

(uo u + uct„J 


11 


12. 


f 1 = ev, 

p = \£ + 2w 


A = ±1, 
a = =F^ 2 


g 8 = e^(A(uc3 u + u^) 
e^u^ + tQs 






+au In u 


A^O, 
v 2 + Xa 
= 1 


X ± = e» ±L (X(ud u 
+vd v ) 

+{u±l)ud v ) 


1 if 

a = 






A^O, 
v 2 + Xa 
= -1 


e vt (X cos t(ud u + u9„) 
+ (z/cost — sin t)ud v ), 
e ut (X sin t(ud u + 
+ (z/sint + cost)ii(9,;) 





The classification of systems of coupled reaction-diffusion equations has been 
ended. 



8 Discussion 

We complete the group classification of systems of reaction-diffusion equations 
started in papers pQ, [2], where systems with a diagonal and square diffusion 
matrix are studied. 

The case of triangular diffusion matrix considered in the present paper ap- 
pears to be rather complicated mainly due to very large number of versions with 
different symmetries. Indeed, we indicate 54 non-equivalent equations (J2J) with 
an invertible triangular diffusion matrix which are presented in Tables 2-6 and 
68 equations with a nilpotent diffusion matrix which are present in Tables 8-10 
and partly in Tables 2, 3. In addition, 14 classes of equations including first 
derivatives w.r.t. variables x v are collected in Table 7. 

In Tables 2, 3 equations defined up to arbitrary functions are presented. In the 
third column of Table 1 we present the main symmetries of the related equations 
(J2J) which do not admit extended symmetries. 

In Table 3 we present equations (J2J) which admit additional equivalence trans- 
formations (AET) (}39|) (indicated in the fourth column) and possible extended 
symmetries indicated if necessary in the third column after the symbol &. 

Tables 4-6 includes the results of classification of equations (J2J) which are 
defined up to arbitrary parameters. 

The items of Tables 2 and 3 marked by the asterisks are related for equations 
((2j) with a ^ and a = as well. 

Tables 7-10 present the results of group classification of equations (j2J) with 
the nilpotent diffusion matrix only. 
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Note that we did not consider linear equations whose group classification is a 
rather trivial problem. 

Among the classified equations there are only seven of them being invariant 
w.r.t. the Galilei transformations. In accordance with Item 2 of Table 3 the 
general form of Galilei-invariant equations ^ is 

u t -aAu = uF 1 , , . 

v t — Au — aAv = uF 2 — vFi 

where F\ and F 2 are arbitrary functions of variable £ = ite«. In Items 3, 5 
of Table 4 and Items 1-4 of Table 6 all functions Fi, F 2 are presented which 
correspond to various possible extensions of the Galilean symmetry. The only 
system of equations (J2J) which admits the extended Galilei group including the 
dilatation and conformal transformations is given by formula (|44J) . 

Thus we end the group classification of systems of coupled reaction-diffusion 
equations (JTJ) started in paper [3] and continued with varying success in 
and PP, The number of non-equivalent equations of this type appears to be 
enormously large (> 300). Nevertheless, using the approach presented in pQ it was 
possible to make an effective classification of pairs of coupled reaction-diffusion 
equations with general diffusion matrix. Moreover, we classify the equations with 
arbitrary number of independent variables. 

We notice that group analysis of react ion- diffusion equations is being inten- 
sively developed in many lines including equations with arbitrary elements de- 
pending on t, x,u,u t ,u x , ■ ■ ■ , refer to [TH] for a survey. Thus our analysis of 
systems of diffusion equations is nothing but a part of general study of diffusion 
models which is currently rather popular. 
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